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We present the results of a calculation of the positive-parity ground-state charmed-baryon spec- 
trum using 2-1-1-1-1 flavors of dynamical quarks. The calculation uses a relativistic heavy-quark 
action for the valence charm quark, clover- Wilson fermions for the valence light and strange quarks, 
and HISQ sea quarks. The spectrum is calculated with a lightest pion mass around 220 MeV, and 
three lattice spacings (a « 0.12 fm, 0.09 fm, and 0.06 fm) are used to extrapolate to the continuum. 
The light-quark mass extrapolation is performed using heavy-hadron chiral perturbation theory up 
to O(m^) and at next-to-leading order in the heavy-quark mass. For the well-measured charmed 
baryons, our results show consistency with the experimental values. For the controversial J = 1/2 
Sec, we obtain the isospin-averaged value Ms^^ = 3595 (39) (20) (6) MeV (the three uncertainties 
are statistics, fitting-window systematic, and systematics from other lattice artifacts, such as lattice 
scale setting and pion-mass determination), which shows a 1.7 ct deviation from the experimen- 
tal value. We predict the yet-to-be-discovered doubly and triply charmed baryons H*,., Qcc, f2*c 
and Qccc to have masses 3648(42)(18)(7) MeV, 3679 (40) (17) (5) MeV, 3765(43)(17)(5) MeV and 
4761 (52) (21) (6) MeV, respectively. 



I. INTRODUCTION 



In recent years, interest in charmed-baryon spectroscopy has resurfaced. This excitement has been partly triggered 
by the first observation of a candidate doubly charmed baryon S+ (3520) by SELEX [T], as well as its isospin partner 
S++(3460) [g. The SELEX Collaboration later confirmed their observation of S+ (3520) 0, but the BABAR [3], 
BELLE [5], and FOCUS [6] experiments have seen no evidence for either state of the isospin doublet (S+ , S++). The 
SELEX evidence for this doublet implies unprecedented dynamics, the most surprising of which is the 60- MeV mass 
difference between the two states. All other previously observed isospin partners have mass differences one order of 
magnitude smaller. The ground-state doubly charmed baryon has been previously studied theoretically via various 
methods, including: the nonrelativistic quark model [7|, the relativistic three-quark model jl], the relativistic quark 
model ,9J, QCD sum rules [10 , heavy-quark effective theory [11], the Feynman-Hellmann theorem [12], and lattice 
quantum chromodynamics (LQCD) [T3HI9] . Overall, theoretical predictions for this state suggest the Sec mass to be 
100-200 MeV higher than that observed by SELEX. 

There remain many undiscovered doubly and triply charmed baryon states. The recently upgraded Beijing Electron- 
Positron Colfider (BEPCII) detector, the Beijing Spectrometer (BES-III), the LHC, and the future GSI project, the 
antiProton ANnihilation at DArmstadt (PANDA) experiment, will help further disentangle the heavy-baryon spectrum 
and resolve puzzles like the one mentioned above. LQCD calculations serve as direct first-principles theoretical input 
for these experiments. 

Currently, LQCD provides the best option for performing reliable calculations of low-energy QCD observables. 
LQCD is a numerical calculation of QCD, which is necessarily performed in a finite discretized and Euclidean space- 
time volume. These approximations introduce an infrared cut-off (the spatial extent L) and an ultraviolet cutoff 
(the lattice spacing a). The latter of these artifacts has been a source of large systematic errors in the heavy-quark 
sector of QCD. For heavy-quark masses satisfying amq <C 1, it is natural to control the discretization errors using 
nonrelativistic QCD (NRQCD) pD]. NRQCD has proven particularly useful when studying physics regarding the 
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bottom quark, but for lattice spacings < 0.12 fm the charm-quark mass is too small to make the NRQCD approx- 
imation justifiable. Alternatively, one can implement relativistic heavy-quark actions |21H25j . where all ©((amg)") 
corrections are systematically removed. 

Several groups have performed lattice charmed-baryon calculations using the quenched approximation |18l 1191 1261 
127] . Although these have served as benchmark calculations of the charmed-baryon sector, the quenched approximation 
is a large source of systematic error that is difficult to estimate. Three previous groups have studied the charmed- 
baryon spectrum using dynamical quarks [T3HI3 HH US] ■ 

Na et al. [ini [13 performed a rather extensive calculation of charm- and bottom-baryon masses at three different 
lattice spacings (a ~ 0.15 fm, 0.12 fm, and 0.09 fm). They used chiral perturbation theory (xPT)-inspired polynomial 
extrapolations of the light-quark masses but refrained from performing a continuum extrapolation of their results. 
From their results for the doubly charmed baryons, one could infer a 40-100 MeV systematic error associated with 
discretization effects. 

Liu et al. |13H15) did a rather nice exploratory calculation over four different pion masses and performed what is 
probably the best (to this day) chiral extrapolation of the J = 1/2+ charmed-baryon spectrum using a relativistic 
heavy quark action for the charm quark. There are a few places where this calculation could be further improved. 
First, the lightest pion used in their calculation was about 290 MeV; with advances in technology, we can get closer 
to the physical point. For baryons with no light degrees of freedom, this is a minor issue, but for isodoublet doubly 
charmed baryons the light-quark mass dependence is nontrivial. Second, they performed all calculations at a single 
coarse lattice spacing, a ~ 0.125 fm, which lies near the upper limit of reliable spacings for studying charm physics. 
In their work, they used power-counting arguments to give estimates of the discretization effects. In particular, in 
the doubly charmed sector, they assigned a rather conservative systematic uncertainty associated with discretization 
effects, 6Mh^^ = —78 MeV. This is by far their largest uncertainty across all states; for example, their result for 
the lightest doubly charmed baryon is M^^^ = 3665(17)(14)t7g. Lastly, they restricted themselves to studying the 
J = 1/2+ sector. The J ~ 1/2+ and J 3/2+ charmed baryons are related by heavy-quark symmetries, which 
results in their chiral extrapolations being coupled. This is particularly relevant when performing a xPT-motivated 
extrapolation of the (S^cS*^) doublet to the physical point. 

The European Twisted-Mass (ETM) Collaboration recently presented results for Ac, Sc, S*, Sec, 2*2, and ilccc, 
using Nf = 2 dynamical sea quarks with a lightest pion mass of about 260 MeV and a relativistic action for the 
valence charm quark [29] . They used ^PT- inspired polynomials for the light-quark mass extrapolation, neglecting 
©(l/mg) corrections and chiral-log contributions. Despite having performed calculations at three lattice spacings 
a G {0.056(1), 0.0666(6), 0.0885(6)} fm, they refrained from extrapolating their results to the continuum. Instead, they 
estimated a discretization error that is incorporated into their systematics. Although historically, the use of Nf = 2 
dynamical sea quarks was a reasonable approximation, this (like full quenching) introduces a source of systematic 
error that can only be quantified when results are directly compared to Nf = 2 + 1 or Nf = 2 + 1 + 1 calculations. 

In order to confidently deal with systematics due to discretization effects, it is necessary to perform calculations 
with highly improved actions, relativistic heavy-quark actions, and multiple lattice spacings in order to extrapolate to 
the continuum. With these goals in mind, we evaluated the positive-parity ground-state charm-baryon spectrum using 
two pion masses (with a lightest around 220 MeV) and three lattice spacings (a « 0.12 fm, 0.09 fm, and 0.06 fm). 
In this work, we made three extensions to our previous preliminary calculation [30]. Firstly, we used an ensemble 
at the super-fine a « 0.06 fm lattice spacing in order to further constrain the continuum extrapolation. Secondly, 
when extrapolating the charmed-baryon masses to the physical ttItt, we used heavy-hadron xPT (HHxPT) [3TH53] 
at next-to-leading order (NLO) in and in the heavy-quark mass expansion, while in our previous work we had 
restricted ourselves to the LO dependence. In order to do this, we extended previous HHxPT results [351 [3S] to 
include ©(l/mg) corrections. Thirdly, we quantified systematics associated with finite-volume effects, scale setting, 
the determination of mjr, 0{m^,a^mTr) corrections to the expressions used to extrapolate to the physical point, and 
the strange-mass tuning. 

This paper is structured as follows. In Sec. |llj we outline the formulation of the lattice calculation, including 
the actions used for the sea, valence light, and valence charm quarks, as well our procedure for setting the scale 
independently, and the construction of our correlation functions. In Sec. |III[ we present the tuning of the charm-quark 
action and show the results for the charmonium spectrum. In this section, we present the results for the mo^— ^n^^ji 
splitting, which is shown to have rather large lattice-spacing dependence, but the result presented is in agreement with 
experiment when extrapolated to the continuum. Section IV outlines our analysis of the charmed-baryon spectrum 
and includes a detailed discussion of the C'(m^,l/mQ) HHxPT expressions for the masses. In this section, the a 
dependence of the charmed-baryon sector is discussed, as well the systematics mentioned at the end of the previous 
paragraph. Finally, in Sec.|V]we give a summary of our results and a comparison of the yet-to-be-discovered masses 
across different models. 
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6.30 0.09845(9) 0.22670(12) 0.09950(53) (23) 0.22464(27)(35) 
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CI 


6.72 0.09444(9) 0.16204(11) 0.09444(38)(9) 0.16086(29)(68) 


48^ 


X 144 


4.5 


330 


1320 



TABLE I: Details of the configurations and propagators used in this work. The subscript "sea" labels the lightest 
sea pseudoscalar masses from the HISQ action [37j|38j, while the subscript "val" labels the valence masses. The sea 
hadron masses have a single uncertainty due to the statistics, while the valence masses include statistical and 



systematic uncertainty due to fitting- window selection as defined in Sec. IIB Additionally, listed are the spatial [L) 
and temporal extents (T) in lattice units, the value of m-^L^ the number of configurations, and the number of 
measurements performed for each ensemble. 



II. LATTICE FORMULATION 
A. Light-Quark Action 

In this work, we used N f = 2 + 1 + 1 gauge configurations that were generated by the MILC Collaboration with the 
highly improved staggered quark (HISQ) [ST-IT action for the sea quarks. The implementation of the HISQ action, 
first proposed by the HPQCD/UKQCD Collaboration ^39lHT] . has been shown to further reduce lattice artifacts as 
compared to the asqtad action 37J. Staggered actions reduce the number of doublers to four "tastes", which are 
reduced to the desired number of true flavors by taking the fourth-root of the fermionic determinant. As a result, 
staggered actions have two sources of discretization errors. The first is due to the discretization of the derivative, 
while the second is associated with taste-exchange interactions in quark-quark scattering. It has been shown that the 
latter type of errors are suppressed at ^ 1% level when the HISQ action is used for light quarks at lattice spacings of 
0.1 fm or less |41j. Furthermore, its suppression of ©((am)^) errors makes the HISQ action a desirable candidate for 
studying charm physics on the lattice [H]. Lastly, despite the HISQ action being significantly more computationally 
expensive than the asqtad action [42], it is still more economical than a non-staggered action. This has allowed the 
MILC Collaboration to recently generate multiple Nf = 2 + 1 + 1 HISQ ensembles, with a range of lattice spacings 
a E [0.045, 0.15] fm and three light-quark (up, down) masses corresponding to m^r G {140, 220, 310} MeV. This variety 
of ensembles allows for clean extrapolations to the physical pion mass and the continuum limit. 

Hypercubic blocking .43i is implemented on the gauge configurations in order to further reduce the ultraviolet noise 
from the gauge field. For the valence light (up, down and strange) fermions a tree-level tadpole-improved clover- Wilson 
action is used[^ since the construction of baryon operators with staggered fermions is rather complicated. However, 
for the coarser and lighter pion mass ensembles (such as 140-MeV pion mass at 0.12 fm), one runs into the problem 
of exceptional configurations where the clover-Dirac operator has near-zero modes |47| . Thus, in this work, we were 
limited to heavier light-quark masses which correspond to m^r G {220,310} MeV with lattice spacings of around 0.06, 
0.09 and 0.12 fm. 

Because the actions used for the sea and valence quarks differ, the calculation presented here is partially quenched 
and violates unitarity. In order to restore unitarity, it is necessary to match the valence- and sea-quark masses, as 
well as to extrapolate the results to the continuum. Due to the four-fold degeneracy of the staggered action, in the 
continuum limit it has an SU{8)l <E) SU{8)ji (g) U{1)y chiral symmetry. In this limit, each pion obtains 15 degenerate 
partners. A finite lattice spacing breaks this symmetry and lifts the degeneracy f48j. Therefore, there is an ambiguity 
when tuning the valence-quark mass to the sea-quark mass. We chose to simultaneously tune the light- and strange- 
quark masses to assure that the valence pion and kaon masses match those of the lightest Goldstone Kogut-Susskind 
sea pion and kaon masses, as shown in Table [H The Goldst one Kogut-Susskind sea pion is the lightest pion, the only 
one that becomes massless in the chiral limit for a nonzero lattice spacing. Ideally, one would want to perform all 
calculations at a range of light, strange and charm masses and simultaneously extrapolate all masses to their physical 
values. Due to limited in computational resources, we performed calculations at a single strange quark mass, but 



The light clover propagators were generated and provided by the PNDME Collaboration I44H46| . 
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as will be discussed in Sec. II D our determination of rriK at the continuum and physical rUj^ is in agreement with 



experiment. This gives us confidence that the strange-quark mass is tuned properly. 

B. Correlation Functions and Fitting Method 

Before discussing the tuning of the charm-quark action, let us explain how we constructed our correlation functions 
id extrac 
functions 



(i) 

and extracted hadronic masses. For a given interpolating hadron operator, O)^ , we construct the two-point correlation 



X 

where the superscripts i and j label the smearing type of the annihilation and creation operator, respectively, {xo,to} 
labels source location, and {x, t} the sink location. In order to reduce statistical noise, the two-point functions are 
averaged over four source locations for each gauge configuration. 

Both the baryonic and mesonic correlation functions are calculated with gauge-invariant Gaussian-smeared (S) 
sources and point (P) sinks. For the mesons, we use the generalized Prony-Matrix (PM) method [IH] over the 
smeared-smeared (SS) and smeared-point (SP) correlation functions. The PM method uses the fact that each choice 
of smearing parameters corresponds to a particular linear combination of the exponentiated masses {rrij) and the 

corresponding overlap factors (Aj), c'j^\t) = AQ*''e~™°* -I- A^^-'e"™^* -I- • • • . By computing correlation functions with 
two sets of smearing parameters, we can determine the two lowest energy states that have overlap with the interpolating 
operator used by solving the eigenvalue equation 

MyH{t + l)-VyHit) = (2) 

where yjj{t) = (cj^^^(t), C^^\t)). One solution to this equation is given by gS] 



M = 



T+tw 



J2 yHit + i)yJiit) 



V = 



T + tw 



E yHit)vl{t) 



(3) 



where the window size t\Y must be > 1 in order to ensure the matrices within the brackets are invertible. For each 
hadron, r is chosen in order to maximize the plateau of the ground state. The statistical uncertainties of the extracted 
hadron masses are evaluated using the jackknife method. 

We test the PM method for a subset of the baryonic masses and compare the results with those extracted from a 
single-exponential and double-exponential fits to the SP correlation function at large Euclidean time. We find these 
to be in agreement within our systematics, with the single-exponential having the smaller uncertainty. As a result, 
we choose to extract all masses from the single-exponential behavior of the SP correlation function. 

For all energies extracted, we determine the statistical uncertainty and a systematic associated with choosing a 
fitting window [ti^tj]. In order to estimate the latter, for all fitting windows that fall within \ti^tf + 2] we calculate 

the energy, x^, and goodness of fit Q{d) (defined as {2'^/'^T{d/2))^^ dxo{XoY^'^~^^~^°^'^)^ which depends on the 
number of degrees of freedom d and is optimally near 1. From this ensemble of energies, we define the systematic as 
the standard deviation of the energies weighted by Q{d). 

C. Charm-Quark Action 

Since the charm-quark mass is too light to justifiably implement a nonrelativistic action for the lattice spacings used 
in our calculation, it is necessary to use a relativistic action. To systematically remove the 0((r7ica)") discretization 
artifacts (where mc is the charm-quark mass), we use the following relativistic heavy-quark action for the valence 
charm quark 



where Qx is the heavy-quark field at the site x, 7^ are the Hermitian gamma matrices that satisfy the Euclidean Clifford 
algebra a^^ = i[7^,7i,]/2, is the first-order lattice derivative, and F^i, is the Yang-Mills field-strength tensor. The 
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a [fm] ama anid amc2 


Al 


0.11926(77)(51) 1.0291(56)(37) 0.901 0.872 


A2 


0.11926(77)(51) 1.0192(31)(21) 0.900 0.853 


Bl 


0.0871(10)(5) 0.7562(81)(52) 0.561 0.536 


B2 


0.0871(10)(5) 0.7463(52)(25) 0.552 0.522 


CI 


0.0578(13)(19) 0.5148(17)(39) 0.319 0.309 



TABLE II: The lattice spacings and fl masses cited include the statistical and systematic uncertainties due to the 
fitting window. The lattice spacings are determined by the chiral extrapolation of the Q mass to the physical value 
of (rriT^/mrtY for each value of Additionally listed are the two bare masses of the valence charm quarks used for 
each ensemble. 

Note that we calculate the f2 mass [arriQ = 0.5007(65) (96)] on 200 configurations for a 0.06 fm, 220- MeV to fix the lattice spacing 
for ensemble CI. 



parameters {mo, j^, cb, ce} must be tuned to assure ©((mco)") terms have been removed. For the coefficients cb and 
ce we use the tree-level tadpole-improved results [TBHT^ [50] cb = v/uq.ce = 1 + ^^/(2^0) with the tadpole factor 

uo defined as uq =^ (1/3) (X]pTr(J7p) ) , where Up is the product of gauge links around the fundamental lattice 
plaquette p. 

The coefficients niQ and v were simultaneously determined nonperturbatively by requiring the ratio mjg/mn = 
{nirj^ + 3mj/^)/{4:m^) to be equal to its experimental value, 1.83429(56), and {rjc, J/ip} to satisfy the correct disper- 
sion relation, Ej^ = + . In constructing the charmonium correlation functions, we used the local interpolating 
operators shown in Table III The dispersion relation was matched using rjc and J /if) energies at the six lowest mo- 
menta: (0,0,0), (1,0,0), (1,1,0), (1,1,1), (2,0,0), (2,1,0) in units of (27r/-L)a~^, and their rotational equivalents. In 
practice, we performed the initial tuning with a subset of 40 gauge configurations (with four sources each). Clearly 
this procedure does not guarantee correct tuning upon analysis of the full ensemble. Therefore, we used two separate 
charm-quark masses and extrapolated to the physical charmonium mass. These two points allowed us to interpolate 
linearly in arric to the physical charm-quark mass defined by mj^/mn = 1.83429(56). The valence charm-quark 
masses used for each ensemble are shown in Table [ll] Figure [T] shows examples of the resulting dispersion relations 
for the rjc and J/'0 with full statistics after extrapolating to the physical charm mass from one of the ensembles, CI, 
and they show that the slopes are consistent with 1. 



D. Lattice-Spacing Determination and Discussion of mff/mn Ratios 

As mentioned earlier, it is necessary to evaluate the spectrum at multiple lattice spacings in order to simultaneously 
restore unitarity and remove discretization errors. With this in mind, we perform the calculation at three lattice 
spacings, a « 0.06 fm, 0.09 fm and 0.12 fm. For the coarse (a « 0.12 fm) and fine (a w 0.09 fm) lattice spacings, we 
use two different light-quark masses corresponding to k, 220,310 MeV; for the super-fine (a « 0.06 fm) ensemble 
we use a single light quark, m^r ~ 310 MeV. We calculate the VL mass on 200 configurations for a sa 0.06 fm and 
ttLt^ k, 220 MeV to fix the lattice spacing for ensemble CI. 

In order to obtain physical masses in the continuum, it is necessary to determine the lattice spacing for the five 
ensembles used. Currently, the most precise determination of lattice spacings for the MILC ensembles is by the 
HPQCD Collaboration [39]; however, their determinations of the lattice spacings for the B2 and CI ensembles 
remain unpublished. For this reason, we perform our own determination. Due to the small m^-dependence of mo 
(at the few-percent level) we choose to set the scale by extrapolating the value of amn across all ensembles with the 
same value of /3 to the physical pion mass. We define the lattice spacing by dividing arn^^^ by the physical f2 mass, 
1672.45(49) MeV. 

In constructing the correlation functions for the Q,, we use (il)* = e'^^"^ {P'^^'^y^ {<li^^''QT) as the interpolating 
operator. The strange-quark annihilation operator is denoted with color index k, F' = C7* are the symmetric spin 
matrices (where C is the charge-conjugation matrix), P+ = (1 -I- 7'*)/2 is the positive-parity projection operator, and 
{P'^'^y^ = 5'^^ — ^7*7'' are the spin-projection operators for spin-3/2 particles. 
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FIG. 1: A sample dispersion relation for the rjc and J/?/'. The six points correspond to energies (and uncertainties) 
for the at the six lowest-momenta states: (0,0,0), (1,0,0), (1,1,0), (1,1,1), (2,0,0), (2,1,0) in units of (27r/i)a"\ 
and their permutations. The red line is the resulting fit to the data using the relativistic dispersion relation 

= rri^j + c^p^ , and the blue band includes the statistical and systematic errors added in quadrature. The 
energies shown are obtained using the full statistics of the CI ensemble and have been extrapolated to the physical 
charm mass. From the fit we obtain the speed of light and its statistical and systematic uncertainties, 
c^^ = 1.0039(28)(9) and cj/y, = 0.9964(35)(5). 
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FIG. 2: xPT and continuum extrapolations of the kaon mass. The red line indicates the fit of the data that has 
been extrapolated to a = 0. The blue band includes the statistical and systematic errors added in quadrature. 



One can determine mo as a function of via SU (3) xPT, but this expression suffers from rather large expansion 
parameters (rriK/l^x^ ''^v/^x) ^^'^ does not always describe lattice baryon masses well. Alternatively, it has been 
proposed that the hyperon masses can be extrapolated using a two-flavor chiral perturbation theory [5T]. With a 
faster convergence than its three-flavor counterpart, the advantages of this approach are clear. The cost is manifested 
in a larger set of unknown coefficients. Using SU{2) xPT for the hyperons, the mass as a function of to^ up to 
0{ml) is ISlj 



J2) 



(47rA) 



J4) 



log- 



(471/,) 
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log' 



o(6) 



(6) 
■7f2 
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where — 130.7(4) MeV is the pion decay constant, and the {do, /^o, 70} are the low-energy coefficients (LECs) of 
the theory. Because at each lattice spacing we have (at most) two ensembles with two corresponding values of m^, we 
are forced to truncate Eq. [Hjat Oijn^) in order to retain a reasonable level of precision. This truncation introduces a 
systematic uncertainty into our calculations which will be accounted for in Sec. |IV B| 

Further details of the ensembles, including our determination of the lattice spacing and the mass are listed in 
Table|ll] The values determined by the MILC Collaboration are a = 0.1211(2) fm for the coarse and a = 0.0884(2) fm 



FIG. 3; Diagrams that contribute to the charmonium correlation functions. In this work we evaluate the 
contribution from connected diagrams (left) and neglect disconnected diagrams (right). The latter are OZI 
suppressed, and previous lattice calculations have determined their contributions to be consistent with zero 

for the fine. The HPQCD Collaboration performed a rather extensive program in which they determined the lattice 
spacing for each ensemble using three different quantities: Upsilon 25'-lS' splitting, the decay constant of the rjs meson, 
and the ri/a ratio We determine a single lattice spacing for each value of /3 and find central values that are 

consistently below both the MILC and HPQCD central values. This difference in the definition of the lattice spacing 
should have no impact on continuum-extrapolated results. 

Table [n] shows that the lattice spacing for the ensemble CI is currently determined at the ^ 4% level of precision. 
For the same reasons discussed above, we choose to determine the physical hadron masses using the m^f/mo ratio. 
As will be shown, the rri/f/mn is determined at the sub-1% level for all ensembles and particles. Due to the removal 
of the C'(m^)-terms in Eq. [sj we proceed to truncate all of our chiral fits at the 0{m^) level of accuracy, and estimate 



a systematic error associated with this approximation (see Sec. IV B). 

Because we are using the strange mass to set the scale, it is important to first test the strange-mass tuning, which 
we do using the kaon mass. For all the pseudoscalar mesons, we use the standard local operators Oh = qflsQf', where 

is the annihilation operator for a quark with flavor / and color index k. As discussed in Ref. |52j . when reducing 
the symmetry of xPT from SU{3) to SU{2), kaons can be represented as a matter field that couples to the SU{2) 
chiral currents. This treatment of the kaons is referred to as KxPT- The advantage of XxPT is that the largest 
expansion parameter is m'^/mj^ < ni\/ {Air fT^Y . Using iiTxPT, the kaon mass as a function of is found to be ^ 

^-K + ^^\'^+0{ml), (6) 



ran ntn \ ^T^Umn J ma 

where is the bare kaon mass, and the a-dependence is parametrized by Ca{iTi^^^a)^ . For the kaon and for all 
other hadrons studied in this work, the continuum-limit mass is recovered by multiplying the ratio at the physical 
point by the physical mass, 1672.45(49) MeV. 

In Table |l] the valence kaon masses are shown for each ensemble. In Fig. [2] we show the values for the kaon mass 
for each ensemble with the corresponding statistical and systematic uncertainties as a function of TOtt/too, as well 
as the chiral extrapolation at the continuum. Figure [2] shows that the lattice-spacing dependence of the kaon is 
rather small, and that the extrapolated value, = 488.7(5.3)(5.3)(5.7) MeV (the three uncertainties are statistics, 

fitting-window systematic, and systematics from scale setting, 0(m^, a^m^) -corrections to the expressions used to 



extrapolate to the physical point, finite volume, and strange-mass tuning as discussed in Sec. IV B), agrees with 
experiment within our systematics. This confirms our strange-mass tuning as well as our scale determination and 
extrapolation procedure using the m_f//mo ratio. 



III. CHARMONIUM SPECTRUM 



In this section, we calculate the charmonium IS splitting and the rest of the charmonium spectrum in the continuum 
limit, and we compare them with experimental and previous dynamical lattice results. We use the ratios of spin 
averages of rjc and J/tp masses to baryon masses to tune the charm-quark mass for each ensemble; thus, the 
splitting between them is not fixed in our calculations. Any deviations from the well-measured experimental values 
give us an estimation of the final systematics. 

In constructing the meson correlation functions, we restrict ourselves to the local interpolating operators shown 
in Table |III[ In order to evaluate the full correlation functions of the charmonium spectrum, we need to perform 
two different types of propagator contractions, as depicted in Fig. |3j connected and disconnected diagrams. Dis- 
connected diagrams would increase the number of propagators needed by approximately two orders of magnitude 
but are suppressed by the OZI rule [SIT. Previous lattice calculations at zero temperature have shown disconnected 
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t t 

FIG. 4: Sample effective-mass plots of the charmonium spectrum from the various ensembles. The errorbar shown 
includes the statistical and systematic uncertainty (from varying the fitted range) added in quadrature. 



Hadron 


'^'^^'Lj J^"^ Interpolator 


{xciY 


'So 0-+ QirM 
'Po 0++ a*QS 

1++ g.SsYQc 

i+- Qh'l'Ql 



TABLE III: Interpolating operators for the charmonium spectrum. Q*^ labels the the charm quark with color index 
k. 



diagrams in the charmed sector are rather noisy, and their contributions to the hyperfine splitting are in the range of 
1-4 MeV and consistent with zero [551 - f55] . Thus, we neglect contributions arising from disconnected diagrams here. 
Figure [4] displays examples of the effective- mass plots after performing the generalized Prony-Matrix method for the 



charmonium sector, and the charmonium masses for each ensemble are shown in Tabic IV in lattice units. 

For every hadron, we calculated the ratio of its mass to the $7 mass, m/z/msi at the two different values of the 
charm-quark mass. After interpolating these to the physical charm-quark mass for each ensemble, we simultaneously 
extrapolated the five values of the hadron masses to the continuum and the physical m^. To perform the light-quark 
mass extrapolation we use the SU{2) xPT expressions, which up to 0{rn^) is linear in m^: 

= — + . f <^cc + 0{ml), (7) 

TOO 47r/^TOf2 V TOO / 

where m^^ is the bare charmonium mass, and the a-dependence is parametrized by Caim^^'^a)'^ . 



Using this procedure, we have verified that our calculations reproduce the experimental low-lying charmed-meson 
spectrum. In Fig. [5] we show our results for the charmonium spectrum (as well as the hyperfine splitting A15 = 
Mj/tp — Mri^) after extrapolating to the physical point. As a result, our errorbars are larger than those of other 
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Hadron 


nic {amH)Ai {amH)A2 iamH)Bi {amH)B2 iamH)ci 


IS 


rrici 1.86213(61)(21) [8-16] 1.85571(32)(5) [9-15] 1.37703(27)(18) [11-18] 1.36696(41)(26) [13-23] 0.93723(18)(11) [12-32] 
mc2 1.83438(47)(24) [8-16] 1.80666(32) (5) [9-15] 1.34397(28)(20) [11-18] 1.32897(41)(15) [13-23] 0.92157(18)(9) [12-32] 


Vc 


mci 1.86213(61)(19) [7-13] 1.85571(32)(4) [16-24] 1.37703(27)(15) [24-36] 1.36696(41) (39) [20-25] 0.93723(18)(14) [18-36] 
mc2 1.83438(47) (27) [7-13] 1.80666(32)(4) [16-24] 1.34397(28)(20) [24-36] 1.32897(41)(13) [20-25] 0.92157(18)(16) [18-36] 


J/^ 


rrici 1.91025(50)(29) [5-17] 1.90212(54)(21) [13-17] 1.41634(78)(19) [26-30] 1.40612(53)(18) [21-26] 0.96470(29) (38) [18-36] 
1.88354(44)(30) [5-17] 1.85446(55)(15) [13-17] 1.38428(81)(29) [26-30] 1.36975(54)(9) [21-26] 0.94955(30) (22) [18-36] 


XcO 


m^i 2.1382(22)(19) [4-9] 2.1264(23)(34) [5-12] 1.5873(29)(27) [7-10] 1.5599(44)(22) [12-23] 1.0619(34)(20) [17-22] 
mc2 2.1126(19)(17) [4-9] 2.0787(23) (28) [5-12] 1.5537(28)(25) [7-10] 1.5209(46)(21) [12-23] 1.0557(19)(20) [17-22] 


Xcl 


mci 2.164(11)(5) [11-16] 2.1574(56)(36) [9-12] 1.6121(26)(12) [3-9] 1.6001(53)(31) [12-23] 1.0966(37)(16) [12-18] 
mc2 2.133(10)(4) [11-16] 2.1104(57)(44) [9-12] 1.5807(26)(16) [3-9] 1.5631(55)(33) [12-23] 1.0814(39)(10) [12-18] 


he 


mci 2.1612(93)(60) [7-10] 2.1573(54)(62) [9-14] 1.6296(59) (45) [10-17] 1.6078(59)(45) [12-23] 1.0904(89)(36) [18-24] 
mc2 2.1373(90)(65) [7-10] 2.1105(55)(35) [9-14] 1.5952(83) (74) [10-17] 1.5709(59)(47) [12-23] 1.0869(53)(20) [18-24] 


Ds 


mci 1.20785(70)(38)[13-23] 1.20348(65) (29) [8-15] 0.89883(46) (40) [13-29] 0.88914(62)(46)[13-23] 0.61196(53)(37)[18-26] 
mc2 1.19203(69) (40) [13-23] 1.17734(64)(26)[8-15] 0.88112(45) (36) [13-29] 0.86803(59)(45)[13-23] 0.60333(52)(28)[18-26] 



TABLE IV: Charmonium and masses in lattice units for the five ensembles and two charm masses. Errors listed 
are statistical and fitting window systematic. The fitting window is given in square brackets. 



calculations. For comparison, we show in Fig. [5]a sample of previous dynamical lattice calculations that have studied 



the charmonium spectrum. By comparing the level of precision of am^ (see Table |II[) and anicc (see Table IV I, one 
can see that it is the uncertainty of am^ that dominates the overall uncertainty of the nice I ratio. 

The works by Bali et al. and Mohler et al. are far more extensive than the small sample that is being represented 
here. Both groups used the variational method over different sources and sinks to not only extract ground-state 
energies but also those of the excited states. Mohler et al. evaluated the spectrum for the {cc, cs, cl} systems for a 
range of six pion masses ranging from 702 MeV to 156 MeV at a single lattice spacing, a 0.09 fm. On the other 
hand, Bali et al. evaluated the {cc} spectrum, including disconnected diagrams, at three lattice spacings but did not 
provide a continuum-extrapolated result for the spectrum or an estimate of the discretization error. 

The conclusion of Fig. [5] is evident: these non-continuum results come with a large systematic error due to nonzero 
lattice spacing. This error decreases with lattice spacing, but from Fig. [5] it is clear that in order to reproduce 
the physical spectrum, it is necessary to extrapolate masses to the continuum. For example, in the upper figure in 
Fig. [5] we see that despite our masses having the largest uncertainties, ours are the only results that are consistently 
in agreement with experiment. We conclude that previous calculations that do not extrapolate their results to the 
continuum have underestimated their systematic errors. 

When tuning the charm mass to the spin-averaged mass, IS*, the most natural quantity to study is the hyperfine 
splitting A15. As a result, this splitting has received a great deal of attention in the community. One surprising 
feature is that for a finite lattice spacing, A15 is underestimated ^TJ [57]. In our calculations we find the value of 
Ais agrees with experiment only after extrapolating to the continuum. This is consistent with the findings of the 
HPQCD/UKQCD Collaboration [57j and Fermilab Lattice and MILC Collaborations [59], as shown in the lower part 
of Fig. [5] Therefore, it cannot be overstated that charmed quantities need to be evaluated at multiple lattice spacings 
to properly quantify the systematics. 

In order to further test the strange- and charm-mass tuning, we evaluated the A^c = — ?Ti^^/2 splitting. This 
is the binding-energy difference between the heavy-light and heavy-heavy systems; there is no reliable analytical 
procedure for calculating this quantity. Since the strange-charm meson, Dg, has no light degrees of freedom, up to 
0{m^) its mass is linear in m^, therefore the A^j, splitting can be extrapolated using 

^...-'^U(9(.4). (8) 



where A^^, denotes the bare splitting, and we extrapolate to the continuum using Ca{m^^^ aY . 

In Table |IV[ the Ds and 77c meson masses are shown for each ensemble. Figure |6] shows the values for the A^c 
splitting after continuum extrapolation, along with their corresponding statistical and systematic uncertainties (see 
Sec. IV B). Figure [6] shows that the a-dependence of A^c is sizable; in fact, continuum extrapolation is necessary in 



order to reproduce the physical value. In performing the continuum extrapolation of Age, we find the a-dependent 
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a Bali et al. (N^2, a=0.09 fm) 

* Briceno et al. (Nf=2+1+1, a=0 fm) 
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FIG. 5: Our determination of the low-lying charmonium spectrum after extrapolating to the physical point, labeled 
as "Briceno et al.", as well as a survey of previous unquenched lattice calculations [HI [55j [57H59] . Mohler et al. 
determined the spitting between {xcOjXcI, ^c} and IS* 57J; in order to compare their results with ours, we have set 
15' to its physical value, while leaving their hyperfinc splitting unchanged. The statistical uncertainty is shown as a 
thick inner error bar, while the statistical and systematic uncertainties (if estimated in the paper) added in 
quadrature are shown as a larger thin error bar. Our systematic uncertainties include errors originating from the 
fitting window, scale setting, pion mass determination, finite-volume effects, ©(m^, a^m7r)-corrections to the 
expressions used to extrapolate to the physical point, and the strange mass tuning (as discussed in Sec. IV B). The 
orange line/band indicates the experimentally measured masses or splittings with their corresponding 
uncertainties 1601. 



LEG to be Ca = —0.0088(46). Since our determination of the cc and cs spectrum is in agreement with experiment, 
we believe that the estimates of the systematics in Sec. IV B accurately reflect the sources of systematic error of the 
calculation presented in this paper. 



IV. CHARMED-BARYON SPECTRUM 



With confidence that our tuning reproduces the low-lying cc, cs. Is spectrum within our systematics, we proceed to 
evaluate the positive-parity charmed-baryon spectrum. Heavy-quark symmetry dictates that the quantum numbers 
of the light degrees of freedom of any heavy-light system are conserved. One can identify approximately degenerate 
multiplets by these quantum numbers. For singly charmed baryons, the light degrees of freedom can have total spin 
equal to zero or one. Under SU{3)v chiral symmetry, the spin-singlet multiplet transforms as a 3 irrep. The spin 



11 




FIG. 6: xPT and continuum extrapolations of the Asc = 1^0^— 1^x1^/2 splitting. The red line indicates the fit of the 
data that has been extrapolated to a = 0. The blue band include the statistical and systematic errors added in 
quadrature. 




t t 



FIG. 7: Sample effective-mass plots from the various ensembles of the charmed-baryon sector. The errorbar shown 
includes the statistical and systematic uncertainty (from varying the fitted range) added in quadrature. 



triplet is a 6 irrep when the total angular momentum is J = 1/2 and a 6* irrep when the total angular momentum 
is J = 3/2. In the heavy-quark limit, these are degenerate. The doubly charmed baryons form a 3 irrep when 
the total angular momentum is J = 1/2 and a 3* irrep when the total angular momentum is J = 3/2. The triply 
charmed baryons are singlets under SU{3)v- This algebra was manifested by the interpolating operators used in this 
calculation, as shown in Table III [Hi] . Figure [t] displays examples of the effective-mass plots for various correlation 
functions. Table |VI| lists the baryon masses in lattice units for each charm-quark mass and ensemble along with the 
statistical and fit ting- window systematic uncertainties and the chosen fitting window. 
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= e*=""P+(P|/')'^Q,*(g,!fr^gr), 
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(fie)' =e'''"^P+{P^^^] 
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(He,)' =e'=""P+(Py'; 


'^gS (Of r^Or) 




(n^,)' -- 




(t7cc)* = e'=""P+(Pl/' 











TABLE V: The interpolating operators for the positive-parity baryons i61j . qu,d,s respectively denote the up-, down- 
and strange-quark annihilation operators, Qc denotes the charm-quark operator, {fc,Z,m} are color indices, while 

denote polarization indices. (r"^,r*) are the antisymmetric and symmetric spin matrices (C75,C7'), where C 
is the charge-conjugation matrix. In order to have the best possible overlap with the state of interest, we have used 
the spin projection operators {P'^^'^Y^ 
projection operator P'^ = (1 -|- 7*)/2. 



iy7^' and {P^^y^ 



_ [py^yi = ^7*7^, and the positive-parity 



FIG. 8: Two of the self-energy diagrams contributing to the masses of a singly charmed baryon in the 6 irrep. The 
first depicts contributions arising from loops containing a pion and a member of the 6 irrep, while the second 
correspond to loops containing a pion and a member of the 6* irrep. There are similar self-energy diagrams for 
baryons in the 6* irrep. 



A. Chiral and Continuum Extrapolation 



As discussed in Sec. II C the ratios of each charmed-hadron mass to the mass are interpolated to the physical 
charm mass, defined by mjg/mn = 1.83429(56). After this is done for each ensemble, it is necessary to extrapolate 
the ratios to the physical light-quark mass and continuum. Due to the rather large expansion parameter of SU{S) 
xPT and poorer convergence rate, we use SU{2) HHxPT to extrapolate the baryon masses to the physical pion mass. 
Previous HHxPT calculations of the singly charmed-baryon masses used the static limit, ttlq — 7> oo [351 ES]- At 
O^I/uiq) new operators are introduced that explicitly break the 6-6* degeneracy resulting in three independent 
bare sphttings {A3 g, Ag g. , A6,6*}- We extend previous work to include the ©(l/mg) corrections for the {Ac, Sc, S*} 
and {Sc, S^, S*} multiplets by evaluating the contribution arising from the two self-energy diagrams depicted in Fig.js] 

First, consider the {Ac, Ec, E*} multiplet. Up to ©(to^), the dependence of the ratio of the particle masses to 
mn can be written as 



ma 



m 



Ac , (^K'^l 65! /I A ^ , 2 X-/ A ^ 



^ _l ^ — ^ 

+Al^ CTv ml 



[■iTrUjmn 3(47r/^)^mo (47r/^)^mo \9 9 

7r~F~\ q/. ^/^2 -^(™T, -Aa,s; , + r\ f l7•^("^7r, 0, fi) + ^ -/^(m^, -As, 

(47r/^)mf2 3(47r/^)^mj7 (47r/^)^mj7 \9 9 



where ajj — {<th — mnPjjcrn/m^^ wP and A" label the bare masses and splittings, and y's and ti's are the LECs of the 
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Hadron 


rUc {amH)Ai {amH)A2 {amH)Bi {amH)B2 {amH)ci 


A, 


fUci 1.4561(42)(70)[8-16] 1.4228(77)(73)[9-15] 1.0808(42)(33)[11-15] 1.0328(102)(79)[16-25] 0.7339(56)(15)[19-22] 
mc2 1.4401(42)(70)[8-16] 1.3976(76)(69)[9-15] 1.0643(41)(35)[11-15] 1.0136(98) (62) [16-25] 0.7258(56) (15) [19-22] 


He 


mci 1.5333(24)(28)[8-15] 1.5120(31)(20)[8-15] 1.1438(37)(21)[14-18] 1.1115(37)(47)[14-27] 0.7747(48)(10)[26-29] 
mc2 1.5174(24)(27)[8-15] 1.4871(31)(21)[8-15] 1.1274(37)(19)[14-18] 1.0922(33) (20) [14-27] 0.7665(48) (10) [26-29] 




mci 1.5521(40)(30)[8-11] 1.5286(50)(54)[8-16] 1.1703(43)(25)[11-16] 1.1351(80)(78)[13-23] 0.7968(32) (54) [13-22] 
1.5359(40) (30) [8-11] 1.5028(50)(51)[8-16] 1.1530(43)(29)[11-16] 1.1134(74)(52)[13-23] 0.7883(32) (54) [13-22] 


K 


m^i 1.6178(43)(48)[7-11] 1.5760(91)(44)[9-15] 1.1979(83)(53)[13-19] 1.1731(105)(167)[13-20] 0.8055(83) (29) [19-24] 
mc2 1.6020(43)(50)[7-11] 1.5516(91)(42)[9-15] 1.1812(82)(55)[13-19] 1.1569(97)(76)[13-20] 0.7975(83) (29) [19-24] 


' — 'c 


mci 1.5878(60)(78)[12-23] 1.5820(55)(54)[11-18] 1.1925(51)(14)[16-22] 1.1682(49)(34)[15-21] 0.8089(23) (22) [12-23] 
mc2 1.5717(60)(86)[12-23] 1.5564(55)(51)[11-18] 1.1753(50)(15)[16-22] 1.1471(44)(16)[15-21] 0.8005(22) (23) [12-23] 


' — 'C 


mci 1.662(3)(14)[8-18] 1.6388(58)(41)[10-14] 1.2314(65)(41)[15-21] 1.2060(54)(48)[14-21] 0.8328(54) (17) [20-24] 
mc2 1.646(3)(14)[8-18] 1.6142(57)(39)[10-14] 1.2157(64)(39)[15-21] 1.1896(51)(9)[14-21] 0.8248(53) (17) [20-24] 




fUci 1.6487(69)(16)[16-24] 1.6393(22)(24)[8-14] 1.2280(45)(17)[19-23] 1.2129(28)(3)[15-19] 0.8341(25)(25)[18-24] 
mc2 1.6322(69)(16)[16-24] 1.6138(22)(24)[8-14] 1.2112(45)(16)[19-23] 1.1919(25)(3)[15-19] 0.8262(24) (24) [18-24] 




m^i 1.6960(38)(52)[ll-20] 1.6882(27)(29)[8-14] 1.2567(64)(34)[19-26] 1.2493(32)(17)[14-19] 0.8567(24) (24) [15-30] 
mc2 1.6805(38)(52)[ll-20] 1.6638(27)(28)[8-14] 1.2408(64)(29)[19-26] 1.2313(29)(7)[14-19] 0.8489(23) (25) [15-30] 


cc 


mci 2.2349(33)(42)[ll-25] 2.2194(67)(61)[15-22] 1.6628(21)(13)[6-16] 1.6413(46)(17)[17-25] 1.1298(25)(12)[19-29] 
mc2 2.2037(33)(39)[ll-25] 2.1701(66)(56)[15-22] 1.6394(48) (50) [6-16] 1.6070(39)(21)[17-25] 1.1139(25)(12)[19-29] 


' — 'CC 


rrici 2.3053(26)(27)[8-16] 2.2455(115)(72)[15-19] 1.6381(55)(47)[18-26] 1.6801(66)(37)[17-22] 1.1570(91)(32)[32-41] 
mc2 2.2744(25)(27)[8-16] 2.1970(114)(73)[15-19] 1.6808(29)(44)[18-26] 1.6459(56)(27)[17-22] 1.1416(91)(34)[32-41] 




mci 2.2893(28)(9)[17-26] 2.2739(22) (12) [15-27] 1.7008(18)(2)[18-26] 1.6786(33) (14) [24-28] 1.1562(14)(4)[19-29] 
mc2 2.2580(28) (10) [17-26] 2.2247(21)(12)[15-27] 1.6677(18) (3) [18-26] 1.6417(28)(6)[24-28] 1.1403(14)(4)[19-29] 




m^i 2.3385(66)(29)[11-19] 2.3178(31)(19)[15-23] 1.7331(43)(10)[23-29] 1.7180(38)(23)[20-26] 1.1796(21)(6)[26-31] 
mc2 2.3078(66)(29)[11-19] 2.2694(31)(19)[15-23] 1.7001(43)(9)[23-29] 1.6799(35)(16)[20-26] 1.1641(21)(6)[26-31] 




mci 2.9621(16)(9)[16-24] 2.9466(15) (17) [16-24] 2.1953(15)(7)[32-39] 2.1788(18)(2)[21-28] 1.4921(22)(8)[38-43] 
mc2 2.9161(16)(8)[16-24] 2.8753(15) (17) [16-24] 2.1472(16)(8)[32-39] 2.1239(17)(2)[21-28] 1.4690(23) (4) [38-43] 



TABLE VI: Charmed-baryon masses for the five ensembles in lattice units, statistical and fitting window systematic 
uncertainties, and fitting windows. 



theory. The chiral function F is defined as 



Hm. A, ,) . (A^ - + In (^^^Ta^^T^ j " ^^7^) ^ ^ ) ' ^''^ 

with J-{m, 0, jj,) — 7rm^. From Eq.jo] in the static limit we reproduce the previous results |35| , l36 j . For the extrapolation 
to the continuum limit, we consider the lattice-spacing dependence of 6r{a) = CairriQ^^^a)^ for each baryon within the 
same multiplet to have the same behavior. 

In order to stabilize our fits, we evaluate the splittings {Aa^Sc j A^^s. , Aa^e; } for each ensemble and extrapolate 
them to the physical pion mass with the assumption that their lattice-spacing dependence is suppressed. The resulting 



splittings serve as input to the chiral function in Eq. 11 In addition, when minimizing we require the axial couplings 
to be real, > 0. This requirement assures that the HHxPT Lagrangian is Hermitian, and it reduces the parameter 
space of the minimization routine, thereby resulting in smaller uncertainties while leaving the mean values of the 
extrapolated masses unchanged. The scale fi is set to 700 MeV; we do not observe a difference in the results when ji 
is varied among {600 MeV, 700 MeV, 800 MeV}. Using the physical value of m^^y fP^y^" = 12.796(37), we find the 
LECs shown in Table VII In Fig. |9] we display our fits at the continuum (a = 0) along with the value of rriH /nifi for 



each ensemble as a function of mj^/rriQ. From Fig. [9] one sees all masses are within 1.1 cr of the experimental values. 



From Table VII it is evident that only the leading-order term in the chiral expression is determined well. 



14 




FIG. 9: NLO HHxPT and continuum simultaneous extrapolations of {A^, Sc, S*} masses. The red line depicts the 
fit of the data that has been extrapolated to a = 0. The blue band includes the statistical and systematic errors 
added in quadrature. 



mljm°n K^^Jml K^^^Jml 



d.o.f. Q 



1.352(33) 0.112(30) 0.162(72) 1.3(1.7) 1.2(5.2) 1(15) 0.2(4.9) 0(16) 0.0042(71) 6.4 6 0.4 



TABLE VII: Results of SU{2) HHxPT LECs from fits of the {A^, S^, S*} muhiplet masses, 
degrees of freedom, and the goodness of the fit Q{d) (as defined in Sec. IIB ). 



the number of 



Next consider the multiplet {.i^c, S^, s*}: 



m" 



9l 



(47r/,)mn 2{AnUYmn \i 
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J"(m^, AH^H^,/i) + -J"(to^, Ah,h;,m) 



^0 



9l 



m 



mg 



n 



9l 



(47r/^)mo 2 (47r/^)2TOn 



J"(m^,-AH^H',Ai) 



J"(m^,-AH,B;,M) + ^ 



(4./J^mo ^9-^(-^'0,^) + ^.F(m.,AH,H 
8 (4.^]^ 0' + iHm., -A 



The values obtained for the LECs are shown in Table |VIII[ and Fig. [TO] displays our fits at the continuum. From 
Fig. [To] it is clear that the extrapolated masses are within 1.1 cr from the experimental values. Furthermore, we are 
not able to resolve any lattice-spacing dependence for this multiplet, and the chiral extrapolation is close to a constant 
for and S*. 
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FIG. 10: NLO HHxPT and continuum extrapolations of {5^, 5'^, 5*} masses. The red line depicts the fit of the data 
that has been extrapolated to a = 0. The blue band includes the statistical and systematic errors added in 
quadrature. 



mljml Al^^./ml A|^s;/mg^ ^5. o^i A gj c^l^t'^ d.o.f. Q 

1.477(45) 0.054(63) 0.11(16) 0.73(60) 0.1(6.7) -0.4(5.1) 3.0(7.1) 0.0(6.4) 0.006(10) 5.2 6 0.5 

TABLE VIII: Results of SU{2) HHxPT LECs from fits of the {Sc,S;;,S*} masses. 



For the multiplet {Scc,S*^} we use the previously determined expressions |63| to perform the chiral extrapolation 



mE 



mi 



TOO 



(12) 



— + 



as* 



(4^/0 " 



mlmn 



5 4 

-J"(m^,0,/i) + -J'(to^,-Ah^^h.^,^) 



+ 0{mt). 



The results for the LECs are shown in Table |IX[ and Fig. [TT] displays our fits at the continuum. It is remarkable in 
Fig. 11 that the m^-dependence of ms^^ is surprisingly small compared to that of ms* . From Fig. 11 one can also 



observe that our value of ms^^ is about 1.7 cr above the experimentally observed mass. 



0/OaO/O - - 2 / phys 


d.o.f. Q 


2.147(35) 0.025(24) -0.00002(55) 0.00057(60) 0.00008(52) 0.013(19) 


6.3 4 0.2 



TABLE IX: Results of SU{2) HHxPT LECs from fits of the {S^c, S*J muhiplct masses. 
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FIG. 11: NLO HHxPT and continuum extrapolations of {Sj-c, S*^} masses. The red line depicts the fit of the data 
that has been extrapolated to a = 0. The blue band includes the statistical and systematic errors added in 
quadrature. The dashed line indicates the physical point ni-j^/m^ = 0.083453(25). 



Hadron 






/ phys 


x' 


d.o.f 


Q 




1.612(24) 


-0.49(66) 


-0.005(18) 


0.57 


2 


0.57 




1.670(23) 


-0.78(62) 


-0.005(18) 


1.32 


2 


0.27 




2.206(30) 


-0.27(81) 


0.010(24) 


0.58 


2 


0.56 




2.247(33) 


-0.17(88) 


0.010(26) 


0.81 


2 


0.44 




2.857(38) 


-0.7(1.0) 


0.019(29) 


1.16 


2 


0.31 



TABLE X: LO SU{2) xPT LECs of isosinglet states n^niMcc^lc: and ft. 



Lastly, the SU{2) HHxPT extrapolation formula for all isosinglet states, ilc,^c,^cc,^cc^ and ilccc, is given by 

= — (T + Ta — T~\ ^ <-^("^7r)- (13) 

In Table |Xj we summarize the fitted LECs of the five isosinglet states. Figure [l2| shows the continuum extrapolation 
of the yet-to-be-observed {ilcc^lcj^ccc} states along with the value of the ratio of their masses to mn for each 
ensemble. 



B. Systematics 

In performing the continuum and chiral extrapolation, we added five systematic errors in addition to the fitting- 
window error. The first of these arises from the uncertainty in determining TOtt and the lattice spacing. We derive this 
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FIG. 12: Chiral and continuum extrapolations of {ilcc, ^^co ^ccc} niasses. The red line depicts the fit of the data 
that has been extrapolated to a = 0. The blue band includes the statistical and systematic errors added in 
quadrature. The dashed line indicates the physical point ni-^/m^ = 0.083453(25). 



uncertainty by simultaneously varying and the lattice spacing within their corresponding uncertainties (shown in 
Tables [T| and [TT| respectively) when extrapolating the masses to the physical point. This gives an ensemble of energies, 
and we obtain a systematic uncertainty from the standard deviation of this ensemble. 

The second uncertainty is due to finite-volume (FV) corrections. The dominant finite-volume effects for baryon 
with light degrees of freedom from the FV counterpart of self-energy diagrams depicted in Fig. [8] and in the p-regime 
these scale like e~™^^ /{mT^L) [M]. More specifically, up to an overall 0(1) constant, they can be written as [M] 



FV,; 



87r/2 m^L\n\ 



(14) 



Note, the overall constant depends on the axial coupling, which we have found to be consistent with zero (see 
Tables VII VIII and IX). For hadrons with no light degrees of freedom, FV effects come in at 0{m'^) in the chiral 
expansion, and therefore are further suppressed by a factor of m.^/ A^, where A-^ ~ 700 MeV is the chiral symmetry- 
breaking scale, 



FV,h 



-L|n|r 



(15) 



In Table |XI| we evaluate both of these FV effects for hadrons with and without light degrees of freedom. 

In performing the chiral and continuum extrapolation we have taken into account terms coming in at 0{a^, m^, m^, l/mg) 
and neglected ©(to^jO^to^) terms. In order to account for 0{m^) corrections we add a systematics of the form [36] 



(4^A)3 



(16) 



which contributes at the MeV level for our ensembles. 
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<5m^^''[MeV] 




Smf^ [MoV] 


5m^^[MeV] 




Al 


1.3 


0.6 


2.1 


4.5 


0.2 


A2 


0.5 


0.2 


1.0 


3.2 


0.7 


Bl 


1.3 


0.6 


2.3 


2.8 


0.6 


B2 


0.3 


0.1 


0.6 


2.0 


2.3 


CI 


1.3 


0.6 


2.1 


0.5 


2.1 



TABLE XI: Shown are estimates for the systematic errors for each ensemble. From left to right column, they are the 
systematic errors due to finite-volume effects for baryons with light degrees of freedom, finite-volume effects for 
hadrons with no light degrees of freedom, the truncation of the xPT extrapolation formulas, corrections in the 
MAEFT expansion, and the sea/valence strange-mass mismatch, respectively. 



In general, quantities obtained using mixed action have discretization errors arising from artifacts of both the sea 
and the valence actions. From mixed-action EFT (MAEFT) we know that leading order (LO) these artifacts can be 
parametrized in terms of two quantities, c^Amix aud a^Agea |651I66) . The latter contributes to the LO a-dependence of 
the valence pion mass, which have been accounted for in our continuum extrapolation. Both of these splittings give rise 
to NLO corrections to the MAEFT extrapolation formulas. In Ref. [67], Orginos and Walker-Loud evaluated a^Ajviix 
for the asqtad improved MILC lattices with a ~ 0.125 fm and found it to be (316(4) MeV)^, which is smaller than 
the corresponding value of a^Agca — (450 MeV)^. To this day, o^Amix has not been determined for the HISQ MILC 
lattices. Assuming this hierarchy (a^Agea > o^Amix) is also true for the HISQ lattices, we can use power-counting 
arguments to estimate the 0{a?"mT^) corrections. 

The values of the a^Agca splittings, which is the mass difference between the Goldstone Kogut-Susskind sea pion and 
the staggered taste-singlet meson, have been determined numerically by the MILC Collaboration for the ensembles 



we are using [371 13H]- From these values we obtain the 6m^^ shown in Table XI Note, the overall 0(1) constants 
present in this correction depend on the axial coupling. 

Furthermore, since we have used the strange mass to set the scale, we need to account for possible mismatch 
between the sea and valence strange-quark masses. We use power-counting arguments to estimate the leading-order 
correction: 

s |(m|.)val - (m?,)goa| , , 

where Jk = 156.1(9) MeV is the kaon decay constant. One can certainly include a similar error for the light-quark 
mismatch, but this would be below our level of precision (0.1 MeV). 

We then add these five sources of systematics for each ensemble and extrapolate them to the physical point, which 
is shown as the third uncertainty of the physical masses in Table |XII[ 



V. DISCUSSION AND CONCLUSION 



In this work we presented the first unquenched continuum determination of the low- lying charmed-baryon spectrum. 
The calculation uses a relativistic heavy-quark action for the valence charm quark, clover- Wilson fermions for the 
valence light and strange quarks, and HISQ sea quarks generated by the MILC Collaboration |n7J[55]. The spectrum 
is calculated with a lightest pion mass around 220 MeV, and three lattice spacings (a ~ 0.12 fm, 0.09 fm, and 0.06 fm) 
are used to extrapolate to the continuum. At each ensemble, we interpolate the charm-quark mass to the physical 
one by matching the charmonium IS spin average through the ratio {m^^^^ + 3m^'^) / {Am^^^) = 1.83429(56); the 
rest of the hadron (composed of charm quarks) ratios mn /rriQ are linearly interpolated in anic to the physical charm 
point. 
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FIG. 13: A survey of previous unquenched lattice calculations [T3l - fT7l l29] . along with the results of this paper 
labeled as "Briceno et al." The "a = 0" label denotes calculations that extrapolated their results to the continuum 
limit. The statistical uncertainty is shown as a thick inner error bar, while the statistical and systematic 
uncertainties added in quadrature are shown as a larger thin outer error bar. Our systematic uncertainties include 
errors originating from the fitting window and scale setting. The experimentally determined masses are shown for 
comparison [60] . 



Hadron 


Latt. Pred. [MeV] 


Exp. [MeV] 


Hadron 


Latt. Pred. [MeV] Exp. [McV] 


Vc 


2995(26)(12)(5) 


2980.3(1.2) 




2481(24)(15)(7) 


2454.02(2) 


J/^ 


3092(27)(13)(6) 


3096.916(11) 


K 


2559(30)(15)(7) 


2518.4(6) 




3397(31)(15)(6) 


3414.75(31) 


' — 'C 


2568(25)(12)(6) 


2575.6(3.1) 


Xcl 


3540(38)(19)(5) 


3510.66(7) 


' — 'c 


2655(26)(6)(7) 


2645.9(6) 


he 


3559(37)(18)(6) 


3525.41(16) 




2681(31)(12)(5) 


2685.2(1.7) 


Ais 


110.9(1. 1)(1.4)(5.3) 


116.6(1.2) 




2764(30)(14)(5) 


2765.9(2.0) 




1960(17)(18)(5) 


1968.45(33) 




3595(39) (20) (7) 


3518.9(9) 


Ds - Vc/2 


468.7(4.8)(5.6)(5.8) 


478.30(69) 


' — cc 


3648(42) (18) (7) 




K+ 


488.7(5.3)(5.3)(5.8) 


493.677(16) 




3679(40) (17) (5) 






2291(37)(22)(7) 


2286.46(14) 


mc 


3765(43) (17) (5) 






2439(29)(25)(7) 


2467.8(6) 


c c c 


4761(52)(21)(6) 





TABLE XIL Results for the charmed hadron spectrum after extrapolating the masses in Tables IV and VI to the 
physical point. The first uncertainty is due to statistics, the second to the fitting window error, and the third 
corresponds to scale setting, finite-volume effects, ©(m^, a^m^)-corrections to the expressions used to extrapolate to 
the physical point, and strange-mass tuning errors added in quadrature (as discussed in Sec. IV B). 
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FIG. 14: Comparison of our results (LQCD) for the masses of the hghtest doubly and triply charmed baryons, with 
the theoretical prediction from other models: QCD sum rules (QCDSR) [13 HH]: the nonrelativistic quark model 
(QM) [7], the relativistic three-quark model (RTQM) [F, the relativistic quark model (RQM) [H], heavy-quark 
effective theory (HQET) [U], and the Feynman-Hellmann theorem (FHT) [12]. 



In order to determine the lattice spacing for the five ensembles, we chose to use the fl mass due to its weak m^- 
dependence. This was done by extrapolating the value of amn over all ensembles with the same value of /3 to the 
physical pion mass. We then obtained the lattice spacing by dividing ama by the physical il mass. The resulting 
values of the lattice spacing are shown in Table [llj 

The main result of this work is the charmed hadron spectrum shown in Table |XII[ which was obtained by ex- 
trapolating measurements from the five ensembles to the physical point defined by m^^^/m^^^ = 0.083453(25) and 
a = (60) . When performing the chiral and continuum extrapolation we use HHxPT up to ©(to^, 1/ 171^0^). The 
three uncertainties of the masses shown correspond to statistics, fitting-window error, and systematics from other 
lattice artifacts, such as lattice-scale setting and pion-mass determination (as discussed in Sec. IV B). 

To test our tuning and extrapolation procedure, we verify that our calculation does in fact reproduce the well known 
low-lying Is, cs, cc spectrum. Since we use the strange quark mass to set the scale, we first determine the kaon mass. 
As is shown in Figs. [2] after extrapolating to the physical point we obtain 771^+ = 488.7(5.3)(5.3)(5.7) MeV, which 
is in perfect agreement with experiment and it displays a minimal lattice-spacing dependence. The remaining results 
for the cs, cc spectrum are shown in Figs. [5]and[6j and it is evident that we do in fact recover the physical spectrum 
in the mesonic sector. Two particularly interesting quantities are the Ds-ric/2 splitting and the A15, both of which 
show significant a-dependence. The fact that we only obtain agreement with experiment after extrapolating to the 
continuum confirms the necessity of performing calculations of the charmed spectrum at multiple lattice spacings. 

In Fig. [13] we display the results for the charmed-baryon spectrum, along with a survey of previous unquenched 
lattice calculations [TSVfrTl and corresponding experimental values for comparison [BH] • Liu et al. evaluated 
the charmed-baryon spectrum for four different pion masses (with lowest ~ 290 MeV) and a single lattice spacing 
a ss 0.125 fm. Na et al. [16l|T7j evaluated the charmed-baryon spectrum at three different lattice spacings (a « 0.15 fm, 
0.12 fm, and 0.09 fm) but have yet to present extrapolated values of the masses as well as an estimate of their systematic 
uncertainties. The European Twisted Mass (ETM) Collaboration recently presented determined the masses of Ac, 
Sc, S*, Sec, Sec, and ficcc, using Nj — 2 dynamical sea quarks with a hghtest pion mass of 260 MeV and three 
lattice spacings a G {0.056(1), 0.0666(6), 0.0885(6)} fm [29]. Despite having performed calculations at three different 
lattice spacings, they did not extrapolate results to the continuum. Instead, they used the information gathered by 
using multiple lattice spacing to quantify the discretization error of their calculations. The use of Nj = 2 dynamical 
quarks, introduces a source of systematic error that is hard to quantify and has not been addressed by the ETM 
Collaboration. That being said. Fig. 13 shows that the masses calculated by the ETM Collaboration are in agreement 
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with our results with the exception of the controversial ^cc, where our result is about 1.6 cr above the value obtained 
by the ETM CoUaboration. 

All previous calculations of the charmed-baryon spectrum have been performed with light-quark masses correspond- 
ing to TO^ > 260 MeV, placing our calculation closest to the physical point. Furthermore, ours is the only determination 
of the continuum charmed-baryon spectrum. Perhaps the most pertinent of the results presented is the Sec mass, 
3595(39) (20) (6) MeV. Unhke all previous calculations, we performed a coupled extrapolation of the {Sec, S*e} doublet 
to the physical point. Although, this led to a tob^c closer to the experimentally observed value in comparison to our 
previous work [3D], our mean value of m-s^^ is still above the mass observed by the SELEX Collaboration by about 
76 MeV [TJ [5] and our combined uncertainty for this particle is 44 MeV. Therefore, despite the fact that we see no 
strong disagreement with the SELEX result, our result does not agree with their experimentally observed mass. This 
is in contrast with the recently published result by the ETM Collaboration, ms^^ = 3513(23)(14) MeV [29], which is 
the only unquenched LQCD calculation to be in agreement with the SELEX Collaboration. 

Therefore, it remains true that the S^ requires further investigation both from the experimental and the theoretical 
communities. In particular, from the experimental side it would be desirable to obtain a clear determination of the 
isospin doublet (S+ , S++) masses. Although, the SELEX Collaboration has confirmed their observation of S+ (3520), 
the BABAR [4 , BELLE [S], and FOCUS [Q experiments observed no evidence for either state of the doublet. From 
the theoretical side, we expect to be able to perform calculations closer to or at the physical pion mass in the near 
future, thereby reducing the contribution from lattice artifacts. In Fig. 14 we compare our results for the masses of 
doubly and triply charmed baryons with predictions from theoretical models. In particular, we show results obtained 
using QCD sum rules (QCDSR) [TUJISH], the nonrelativistic quark model (QM) [7], the relativistic three-quark model 
(RTQM) [5], the relativistic quark model (RQM) [3], heavy-quark effective theory (HQET) [Tl], and the Feynman- 
HcUmann theorem (FHT) [T^. Our result for toh^^ is 3595(39) (20) (6) MeV, and from Fig. 14 we estimate the overall 
theoretical prediction for this mass to be 3550-3650 MeV. These figures can guide experimentalists on the quest 
for the doubly and triply charmed-baryon masses. Finally, we predict the yet-to-be-discovered doubly and triply 
charmed-baryon masses S*^, Occ, n^^, n^cc to be 3648(42)(18)(7) MeV, 3679(40)(17)(5) MeV, 3765(43)(17)(5) MeV 
and 4761 (52) (21) (6) MeV, respectively. 
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